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Abstract. Second-order phase transitions in a non-equilibrium liquid-gas model with reversible mode cou- 
plings, i.e., model H for binary-fluid critical dynamics, are studied using dynamic field theory and the 
renormalization group. The system is driven out of equilibrium either by considering different values for 
the noise strengths in the Langevin equations describing the evolution of the dynamic variables (effectively 
placing these at different temperatures), or more generally by allowing for anisotropic noise strengths, 
i.e., by constraining the dynamics to be at different temperatures in ci||- and dx- dimensional subspaces, 
respectively. In the first, isotropic case, we find one infrared-stable and one unstable renormalization group 
fixed point. At the stable fixed point, detailed balance is dynamically restored, with the two noise strengths 
becoming asymptotically equal. The ensuing critical behavior is that of the standard equilibrium model 
H. At the novel unstable fixed point, the temperature ratio for the dynamic variables is renormalized to 
infinity, resulting in an effective decoupling between the two modes. We compute the critical exponents at 
this new fixed point to one-loop order. For model H with spatially anisotropic noise, we observe a critical 
softening only in the dx- dimensional sector in wave vector space with lower noise temperature. The ensuing 
effective two-temperature model H does not have any stable fixed point in any physical dimension, at least 
to one-loop order. We obtain formal expressions for the novel critical exponents in a double expansion 
about the upper critical dimension dc = 4 — dy and with respect to dy, i.e., about the equilibrium theory. 

PACS. 64.60.Ht Dynamic critical phenomena - 05.70.Ln Non-equilibrium thermodynamics, irreversible 
processes - 64.60.Ak Renormalization-group, fractal, and percolation studies of phase transitions 



1 Introduction 

The theory of equilibrium dynamics for critical phenom- 
ena has identified a set of universality classes which de- 
scribe the long-time, long-wavelength behavior of systems 
in the vicinity of critical points [Q . The situation is in this 
respect similar to that found in the study of static critical 
phenomena. However, whereas in the static case the na- 
ture of the universality classes is determined solely by the 
nature of the interactions that fixes the dimensionality n 
of the order parameter space as well as the effective spatial 
dimension d for its fluctuations, in the case of dynamical 
critical phenomena the corresponding universality classes 
are not simply characterized by the static or dynamic in- 
teractions which exist between the relevant system vari- 
ables, but they also crucially depend on the conservation 
laws which are implemented by the dynamics and on the 
very existence of the detailed balance constraint that guar- 
entees relaxation to thermal equilibrium in the long-time 
limit. 

At the quantitative level, the equilibrium dynamics of 
critical systems is usually formulated in terms of coupled 
non- linear Langevin equations. Such equations typically 
include a dissipative term which involves the Hamilto- 



nian describing the static critical properties and a Gaus- 
sian white noise term, which mimics thermal fluctuations 
and the random forces originating from couplings to fast, 
non-critical non-conserved modes. Furthermore, the equa- 
tions of motion may also include purely reversible mode- 
coupling terms which represent dynamical interactions be- 
tween the relevant slow system variables. In this frame- 
work, the conditions for the existence of detailed balance 
are (i) the Einstein relation between the relaxation con- 
stants and the noise strengths, and (ii) the condition that 
the probability current associated with reversible mode- 
coupling terms be divergence-free in the space of slow dy- 
namic variables 1^. The existence of conservation laws in 
turn fixes the precise form of the relaxation terms ^ . If 
the coupled system of Langevin equations obeys these two 
conditions, it can be shown, e.g., by means of the associ- 
ated Fokker-Planck equation, that the long-time steady 
state of the dynamics is indeed characterized by a Gibb- 
sian probability distribution, with precisely the Hamilto- 
nian that describes the static critical properties of the sys- 
tem. These conditions also insure that the dynamic sus- 
ceptibilities reduce to the static ones in the limit of zero 
frequency, and imply the validity of fiuctuation-dissipation 
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relations between dynamical correlation functions and the 
dissipative parts of the response. 

On a more general footing, one can also consider the 
study of dynamical phenomena which are truly non-equi- 
librium in the sense that they do not possess a steady 
state described by a Gibbs distribution, of which the best- 
known example is perhaps the Kardar-Parisi-Zhang equa- 
tion (for d > 1) Q, which describes the curvature-driven 
epitaxial growth of a surface. Other prominent examples 
are driven diffusive systems 1q], models of driven inter- 
faces and growing surfaces [||7depinning transitions |7|, 
and phase transitions from active to absorbing states |3|, 
e.g., in diffusion-limited chemical reactions. 

A question which lies at the interface between these 
two subjects, i.e., between the study of systems with equi- 
librium and far-from-equilibrium dynamics, concerns the 
asymptotic scaling properties of a critical system origi- 
nally in thermal equilibrium, which is however subjected 
to a perturbation that violates the detailed balance con- 
ditions. This issue is also relevant in an experimental con- 
text, because the maintenance of thermodynamic equi- 
librium in a critical system during an experiment is a 
non-trivial task, as due to critical slowing-down the re- 
laxation times become very long. For the interpretation 
of the experimental results, it might thus be important 
to know whether the dynamical system will eventually be 
driven to a genuine non-equilibrium fixed point which is 
characterized by scaling exponents distinct from the orig- 
inal equilibrium critical behavior. A priori, a perturba- 
tion from thermal equilibrium may either imply a viola- 
tion of condition (i) or of condition (ii) above. This latter 
case includes, for example, driven lattice gases in which 
the terms added to the Langevin equation modelling the 
equilibrium diffusive dynamics stem from a global gradi- 
ent such as an energy or mass current through the sys- 
tem 1^. The violation of condition (i) corresponds, on the 
other hand, to the coupling of the system to a local en- 
ergy gradient. This type of perturbation was the subject of 
study in Refs. |9| Jl0| , |ll| , p^ , where the consequences of the 
violation of the Einstein relations in several equilibrium 
dynamical models, largely describing magnetic systems, 
were investigated within the framework of the dynami- 
cal renormalization group. Specifically, two generic types 
of detailed-balance violations were studied there, namely 

(a) coupling the order parameter and additional conserved 
quantities to heat baths with different temperatures, and 

(b) allowing for spatially anisotropic noise correlations for 
conserved variables. 

In this paper, we extend the above studies [p| Jl0| , [ll| , ^ 
to a model which describes the dynamics of a liquid-gas 
phase transition, or the phase separation of a binary liq- 
uid, at its critical point. This is the so-called model H in 
the terminology of Ref . jl| , which incorporates the inter- 
action of the conserved scalar order parameter, a linear 
combination of the mass and energy densities, namely es- 
sentially the free energy density, with the conserved trans- 
verse mass current vector via a reversible mode-coupling 
term These two fields are sufficient to describe the 
critical dynamics although critical effects can also be seen 



in other modes, such as the sound mode p4| . Further- 
more, a renormalization group analysis was applied to in- 
vestigate non-universal properties and crossover behavior 
in the equilibrium dynamics . The original approaches 
to the equilibrium critical dynamics of model H utilized 
self-consistent mode coupling theory More recently, 
Patashinski performed a linear analysis of the equations of 
motion to study the effect of perturbations which induce 
a non- uniform non-equilibrium stationary state [0. 

We use the response functional formalism to ex- 
press the dynamical equations as a path integral which 
represents a very convenient form to develop the pertur- 
bation expansion and subsequent renormalization group 
analysis. We compute the beta functions of the theory to 
one-loop order in the isotropic case where the two dynam- 
ical variables are coupled to heat baths at different tem- 
peratures and also in the case where we allow for spatially 
anisotropic noise correlations for these conserved fields. In 
the first situation, we obtain one stable and one unstable 
fixed point. The stable fixed point is just the ordinary 
equilibrium fixed point, for which detailed balance is dy- 
namically restored, and in whose vicinity the two noise 
temperatures become equal. On the other hand, at the 
genuinely non-equilibrium, but unstable fixed point, the 
ratio between the temperature of the noise coupled to the 
order parameter and that of the noise coupled to the con- 
served current is renormalized to infinity. In this limiting 
case, the critical exponents of the conserved order param- 
eter are simply those of the decoupled diffusive model B. 
This is due to the existence of a unidirectional random 
heat flow from the order parameter heat bath to the mass 
current thermal reservoir which renders the effect of the 
mode-coupling terms in the dynamics of the order param- 
eter negligible. This situation is quite analogous to that 
found in the earlier study of the non-equilibrium Sasvari- 
Schwabl-Szepfalusy (SSS) model Q, although there the 
order parameter in not conserved, and hence follows model 
A dynamics in the corresponding limiting case. Moreover, 
owing to the fact that a conserved critical field always 
relaxes much slower than a diffusive mode, here we find 
no fixed point that would describe the above temperature 
ratio scaling to zero. 

When allowing for spatially anisotropic noise, we ob- 
serve a softening of the dynamics only in the dj^-dimensio- 
nal sector in wave vector space with lower noise temper- 
ature. We show that the renormalized coupling constants 
diverge as the dimension of the subspace with higher noise 
temperature approaches = 0.838454; thus, at least to 
one-loop order, there exists no finite fixed point for any 
physical value of . We then obtain formal expressions for 
the novel critical exponents in a double expansion about 
the upper critical dimension dc = 4 — d|| and with respect 
to i.e., about the equilibrium theory. Here again the 
results are similar to those we have previously obtained in 
the study of the two-temperature non-equilibrium model 
J pT| |. These conclusions are thus in line with our ear- 
lier studies and with previous observations that whereas 
models with a non-conserved order parameter are quite ro- 
bust against violations of the detailed balance conditions. 
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models with conserved dynamics seem to be extremely 
sensitive to this type of perturbations provided they are 
rendered anisotropic 

The structure of this paper is as follows: In section 
^ we present the coupled non-linear Langevin equations 
which define our two non-equilibrium versions of model H, 
namely the isotropic case where the two dynamic variables 
are at different temperatures, and the situation where we 
allow for anisotropic noise correlations, and we show how 
these models can be mapped to a dynamic field theory. 
The action of these functionals is then separated into a 
Gaussian part and a non-linear contribution which is to 
be treated in a perturbation expansion. In section |^, we 
present the renormalization group analysis of the isotropic 
non-equilibrium model H to one-loop order and discuss its 
behavior near the renormalization group fixed points. In 
section ^, after introducing spatially anisotropic noise cor- 
relations, we derive the renormalization group flow equa- 
tions for the resulting effective two-temperature model 
H, and we compute the values of the critical exponents 
in a formal double expansion about the equilibrium the- 
ory, also to one-loop order. We demonstrate that there 
exists no one-loop renormalizationg group fixed point in 
any physical dimension, however. Finally, in section ^ we 
present our conclusions. 



2 Model equations 

In this section, we briefly outline the basic model equa- 
tions for our isotropic and anisotropic non-equilibrium 
generalizations of model H. Following the equilibrium the- 
ory p3[ |, we consider a second-order phase transition for 
the scalar order parameter ^o(x, t) = e(x, i) + {fi — Ts) 
p(x, t) at a liquid-gas (or binary-fluid) critical point, where 
e(x, t) is the energy density, p(x, t) the mass density and 
/I, T and s are, respectively, the chemical potential, the 
temperature and the entropy at equilibrium (we denote 
unrenormalized quantities by the subscript '0'). This or- 
der parameter is dynamically coupled to the transverse 
mass current Jo (x, t), which satisfies V-Jq = 0. Notice that 

the longitudinal current jj} is related to the mass density 

through the continuity equation dtp + V -Jq = 0. The two 
dynamic fields ipo and Jq yield three hydrodynamic modes, 
namely thermal and viscous diffusion. A fourth hydrody- 
namic mode, i.e., the sound mode, is also present in a fluid, 
but for low momenta its characteristic frequency is much 
higher than the frequencies of the shear and energy modes 
and one may thus disregard it to first approximation. 

Next, the static critical properties of the system are 
described by the following Landau-Ginzburg- Wilson free 
energy functional in d space dimensions. 



d'x ^Vo'(x) 



Uq 



\ [V^o(x)]' 



+ ^Vonx)+2J^(x) , (2.1) 



where tq = (T - T°)/T° denotes the relat ive deviation 
from the mean-field critical temperature T^. Note that the 
quadratic term in Jq simply represents the kinetic energy 
of a stationary mass current. Since Jq itself is a non-critical 
variable, the coefficient in front of this term will be weakly 
temperature-dependent only, and can thus be set constant 
in the vicinity of Tc, and eventually absorbed into jg. 

This effective free energy (2.1) determines the equilib- 
rium probability distribution for the vector order param- 
eter -00 and for the mass current Jq, 



,-ifo[VoJo]/'cBT 



(2.2) 



Note that it follows from (p]2) that at the purely static 
level the order parameter "00 (x) and the mass current Jq (x) 
are completely decoupled; moreover, Jo(x) is a Gaussian 
variable, whose contribution to the free energy can be 
readily factored out of the functional integral. The task is 
then reduced to the computation of two independent criti- 
cal exponents, e.g., "q and by means of perturbation the- 
ory with respect to the static non-linear coupling uq and 
by employing the renormalization group procedure, within 
a systematic expansion in terms of e = 4 — d about the 
static upper critical dimension dc = 4. Here, 77 describes 
the power-law decay of the order parameter correlation 
function at criticality, (-0(x)0(x')) cx l/|x — x'|^~^+'', or, 
equivalently, of the static susceptibility, x(q) cx l/q^"'', 
and the exponent v characterizes the divergence of the 
correlation length as Tc is approached, ^ oc \T — Tc\~^ . 
Notice that fluctuations also shift the true transition tem- 
perature Tc downwards as compared to the mean-field crit- 
ical temperature T°, i.e., rpc = Tc — < 0. Since we 
will need to consider (time-dependent) correlation func- 
tions of the dynamical variables, it is convenient to add 
to the free energy functional a term involving external 
sources h and A, whereupon the (generating) functional 
reads H = Hq - J d'^x [/i(x)0o(x) + Jo(x) ■ A(x)]. 

With this free energy functional H, the coupled non- 
linear Langevin equations defining model H read |l3] 



(2.3) 



AoV2(ro- V2)V;o + 



and 



dt 



T 



T [DoV% - DoV'A 
+5oVV^o(?-o - V2)V;o 



(2.4) 



(2.5) 



goViPoh + C]- (2.6) 



Here, go denotes the strength of the reversible mode-coup- 
ling terms, and 77 and ^ represent fluctuating forces with 
zero mean, (7y(x, i)) = 0, (C"(x, i)) = 0. The couplings Aq 
and Dq are, respectively, the thermal conductivity and 
shear viscosity (in appropriate units). T[. . .] denotes a 
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projection operator which selects the transverse part of 
the vector in brackets; in Fourier space T"^(k) = S"^ — 
k"k^ /k"^. The above Langevin-type equations of motion 
are invariant under Galilean transformations, a symmetry 
which can be utihzed in the renormalization procedure 
(see below and the appendix). 

In order to fully characterize the dynamics, we fur- 
thermore need to specify the correlations of the stochas- 
tic forces. Given that both the order parameter and the 
transverse current are conserved quantities, the strength 
of the random forces has to vanish at zero momentum. 
We assume that these random fields are Gaussian corre- 
lated and we write, in general, the second moment of the 
distributions in Fourier space as: 



(?7(k, = 2Ao(k) <5(k + k') 6{,y + v') (2.7) 



and 



(r (k, z.)C^(k', v')) = 2I?o(k) <5(k + k') 5{v + v') 

k-^kl^ 



X 5' 



fc2 / ' 



(2., 



where the transverse projector in ( |2.8| ) insures that the 
random force is in the transverse directions only. For the 
equilibrium model H, the functions Ao(k) and £'o(k) are 
equal to Ao(k) = Aq ksTk"^, I)o(k) = Dq ksTP, i.e., the 
noise correlators satisfy the Einstein relations and the con- 
servation conditions Ao(0) = -Do(O) = 0. 

As stated above, we will consider two choices for the 
noise correlators which do not satisfy the Einstein con- 
ditions and therefore describe non- equilibrium versions of 
model H. Both these choices can be viewed as generic per- 
turbations from the equilibrium situation. In the first case, 
we take: _ _ 

Ao(k)-AoA:2 (2.9) 

and _ _ 

Do(k) = Z?oA:^ (2.10) 

with Ao/Ao 7^ Dq/Dq. Since the equilibrium model H cor- 
responds to the case in which Aq/Aq = Dq/Dq = kgT, one 
may interpret Eqs. (2.9) and ( 2.10] ) as describing a situa- 
tion in which the two dynamical variables, the order pa- 
rameter and the transverse current, are effectively placed 
in contact with heat baths at different temperatures. 
In the second case, we take more generally: 



and 



Ao(k) 



Do(k) 



dI kl 



D^ kj 



(2.11) 



(2.12) 



This choice of the noise correlators corresponds to a situ- 
ation where two spatial sectors of dimensions d|| and d± 
are placed at different temperatures, i.e., where real space 
isotropy is broken |l^. Without loss of generality, we set 
the effective noise temperature higher in the parallel sub- 
space. 

Once the noise distributi on is specified, one c an r epre- 
sent the Langevin equations ( |2.4|) and (2^), with (2/7) and 



( p.8[ ), as a dynamic field theory, following standard proce- 
dures ||l|,0. This resuhs in a probability distribution for 
the dynamic fields i/'o and Jq, 



^[{^o,jo}] « J mm] J m^io}] x 



(2.13) 



with the statistical weight given by the Janssen-De Do- 
minicis functional J = Jhar + <^rci + >^mc + <^c, which we 
divide into an harmonic part Jhar, which one can integrate 
exactly; the purely relaxational (of static origin) and re- 
versible dynamic non-linear terms Jroi and Jmc, which are 
to be expanded in a power series, giving rise to the pertur- 
bation series in terms of Feynman diagrams; and a term 
Jsc which depends on the external sources. 

The harmonic part, in terms of the original dynamic 

fields ?Ao(k, i^), Jo(k, :/) and the auxiliary fields V'o(k, i^), 
^□(k, i/) (we use, for convenience, the Fourier space repre- 
sentation) reads 



Jhar[{V'o},{V'o},{jo},{jo}] 



(2.14) 



{2TrY 2tt •^o(ki)'(/'o("ki,-z^i)V'o(ki,i^i) 
+ J2 5o(ki) j^(-ki, -^.i)r"'5(ki)J^(ki, i^i) 

a,f3 

- Tpoi-ki, -vi) [-ivi + Xokj (ro -I- kf)] ■(/'o(ki, i^i) 

- ^j^(-ki, -i.i)(~*z.i + DQkl)T^P{k,)/,{kui^i) 

a,l3 

where the projector T°'^{'ki) insures that only the trans- 
verse components of the fields j^, jg contribute to the ac- 
tion. The longitudinal component, as mentioned before, 
does not represent independent fluctuations, and its con- 
tributions (though formally infinite) can be factored out 
of the functional integral. 

The static non-linearity leads, in turn, to the relax- 
ation vertex 

J„i[{^o},{^o}] (2.15) 
XqUq f d'^ki d'^k2 d'^k^ dvi dv2 dv^ 



6 



{2ttY {2t:Y {2ttY 2-k 2it 2t: 



X fc^^o(~k, -v) ?Ao(ki, vi) V'o(k2, 1^2) '0o(k3, 1^3), 

where (k, u) — (ki -|- k2 -I- ka, ^1 + + ^^s)- The purely 
dynamic couplings generate the mode-coupling vertices 



Jmc[{V'o},{V'o},{jo},{Jo}] 

d^fci d'^k2 dvi dv2 
(27rp {2iiY '2tt '2tt 



-Wo 



E 

a,0 



fc? r"'3(k2) V'(-k, -i.)V^(ki, i.i) j'^(k2, 1/2) 

-Uk?kl + k^kf)T'^f\k) 



X ,]^(-k, -i^) V'(ki, ;/i)V'(k2, ^2) 



(2.16) 
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with (k, ly) =^ (ki + k2, i^i + 1^2) here. 

Finally, the source terms give rise to the following con- 
tribution 



^sc[{V'o},{V'o},{jo}>{Jo} 
d^fci dvi 



(2.17) 



Aofcj "(/"(-ki, -vi) /i(ki, j/i) 



J2 Doklf{-k,, -v^) r"''(ki) A'^(ki, y^] 
dfiki cfik2 dvi dv2 



{2itY (2TrY 2tt 2tt 



aui au2 s—^ 

OtT 9-77 2^ 



a, 13 



fc? r"'3(k2) ^(-k, -i.)^(ki, ly,) A^ik2, V2 
k'^ T"^(k)/(-k, -V) /l(ki, :/i)7A(k2, :/2) 



with (k, v) = (ki +k2, J^i + (^2) in the second integral. The 
contribution of the sources is simply included for conve- 
nience. If one wishes to relate the dynamic susceptibilities 
to the Green's functions, one can simply differentiate the 
dynamic functional with respect to /i(x, t), A(x, i) and 
then set these sources equal to zero. The exp ression of 
t he d ynamic functional in Fourier space ( ^.14 ) through 
(2.16), although more cumbersome in notation, is directly 
applicable to the perturbation ex pans ion. 

As usual, the harmonic part (2.14) defines the propa- 
gators of the field theory, while the perturbation expan sion 
i s per formed in terms of the non-linear vertices ( 2.15 ) and 
. Notice that the existence of the reversible forces 
does not show up in dynamic mean-field theory (van 
Hove theory), which i n fie ld-theory language is based on 
the harmonic action ( 2.14 ) only. Wc shall see that while 
the choice ( |2.9| ), ( ^.10 ) for the functions Ao(k) and _ Do(k ) 
yields a perfectly consistent field theory, the choice ( 2.11D , 
( 2.I2D predicts a correction to the critical temperature 
which is anisotropic. One therefore needs to modify the 
theory accordingly in order to obtain a consistent descrip- 
tion of the model. We shall address this issue in section 
^. Here it suffices to say that once we have modified the 
Langevin equations describing the dynamics, we can also 
treat the resulting model by the means described above. 

Once the dynamic functional is obtained, the pertur- 
bation expansion for all possible correlation functions of 
the dynamic and auxiliary fields, as well as for the as- 
sociated vertex functions, is given by the one-particle ir- 
reducible Feynman diagrams. A straightforward scaling 
analysis yields that the upper critical dimension of the 
isotropic model H is dc = 4 for both the relaxational and 
the mode-coupling vertices. For the anisotropic, effective 
two-temperature model H, however, the upper critical di- 
mension is reduced to dc = 4 — d\\ , as will be seen be- 
low. For d < dc, the perturbation theory will be infrared- 
singular, and non-trivial critical exponents ensue, while 
for d > dc the perturbation theory contains ultraviolet 
divergences. In order to renormalize the field theory in 
the ultraviolet, it suffices to render all the non- vanishing 
two-, three-, and four-point functions finite by introduc- 
ing multiplicative renormalization constants, following an 



additive renormalization corresponding to a fluctuation- 
induced shift tqc of the critical temperature. This is achie- 
ved by demanding the renormalized vertex functions, or 
appropriate momentum and frequency derivatives thereof, 
to be finite when the fluctuation integrals are taken at a 
conveniently chosen normalization point, well outside the 
singularities of the infrared regime. 

We shall employ the dimensional regularization scheme 
in order to compute the emerging momentum integrals, 
and choose the renormalized mass r = 1 as our normaliza- 
tion point, or, sufficient to one-loop order, tq = ro — roc = 
/U^. Notice that /i defines an intrinsic momentum scale of 
the renormalized theory. From the renormalization con- 
stants {Z factors) that render the field theory finite in the 
ultraviolet (UV), one may then derive the renormalization 
group (RG) flow functions which enter the Gell-Mann- 
Low equation. This partial differential equation describes 
how correlation functions change under scale transforma- 
tions. In the vicinity of an RG fixed point, the theory 
becomes scale-invariant, and the information previously 
gained about the UV behavior can thus be employed to 
access the physically interesting power laws governing the 
infrared (IR) regime at the critical point {t ck T — Tc —> 0) 
for long wavelengths (wave vector q — > 0) and low frequen- 
cies {cj ^ 0). 



3 Renormalization group analysis of the 
isotropic model H 

3.1 Vertex and response function renormalization 

The UV-divergent two-, three-, and four-point vertex func- 
tions or their derivatives that require multiplicative renor- 
malization are S^r^ (q, w) , 8^2 (q, uj) together with 

^i^^oi^i,^'^'^') ^'^'^ ^i^^o^^'^^^'^y^ functions involv- 
ing J, J, a„rQyj(q,w), a,2rg~ (q,w) and dg2r^jj{ci,u}); 

the three-point vertices '9g°^oi/J^i/'J° ^""^ ^ ^' ^' ^' 

and 9(5»p2+po^2)r~»^^^(-q - p, -w; q, w; p, 0); and finally 
the relaxation vertex d„2r~ , , i—a 



'3'3'3'3'3'3''' 

On the other hand, we have four fluctuating fields (j^, Jq , 
ipo, Vo) and the seven parameters Dq, Dq, Aq, Aq, tq, go 
and uq available; this leaves us at liberty to choose one of 
the renormalization constants in a convenient manner. In 
addition to these parameters and since detailed balance 
does not hold for this non-equilibrium model, one also 
needs to consider the renormalization of the dynamic sus- 
ceptibility of the order parameter, in order to determine 
the susceptibility exponent 77. 

Starting with the two-point functions r^~^{ci,Lu) and 

-Tpjj (q, Lu) for the conserved order parameter and trans- 
verse current fluctuations, respectively, we immediately 
note that as a consequence of the momentum dependence 
of the mode-coupling vertices, one has for these two vertex 
functions 

a|j)^o.~.(Tj)(^ = «'-)-i (3.1) 
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to all orders in perturbation theory. Upon defining renor- 
malized fields according to 



J 



Jo 



(3.2) 
(3.3) 



which imply that = (Z-Z^)-^/^ T^-^ and that 
(2~Zj)^^/^r'pjj, we thus obtain the exact relations 



Z^Z,p = 1 , 
At this point wc utilize our fi'ccdom of choice 



(3.4) 
(3.5) 



to set 



(3.6) 



The multiplicative rcnormalization factors for the cou- 
pling constants are defined through 



A = 


Z\\q , 


(3.7) 


A = 


ZjXo , 


(3.8) 


D = 


Zd Do fr"^ , 


(3.9) 


D = 


Z^Do^i-^ , 


(3.10) 


T = 


Zt To with To = To - roc , 


(3.11) 


U = 


ZuUoAdu'^^'^ , 


(3.12) 


9 = 


71/2^ ^1/2 d/2-3 


(3.13) 



where the ^ factors have been introduced to render the 
renormalized couplings dimensionless. The geometric nu- 
merical factor Ad = r(3 - d/2)/2'^- V^^^ has been ab- 
sorbed into u and g since it appears in the explicit inte- 
grals of the subsequent perturbational analysis. Finally, 
we still have at our disposal a symmetry of the theory 
which is valid even in a non-equilibrium situation, namely 
Galilean invariance (see Appendix ^) . This symmetry im- 
poses the exact condition that 

ZgZj=l. (3.14) 

From this relation and Eq. (BM it follows that the dy- 



namical vertex does not renormalize, i.e. Zg = 1. 

In order to discuss the RG flow diagram, it is con- 
venient to introduce the following rescaled coupling con- 
stants 



■"0 = T- "0 
Ao 



Ao 
Ao ffo 



Ao AqI^o 



(3.15) 
(3.16) 



and the ratio between the noise temperature of the order 
parameter and the noise temperature of the mass-energy 
current, respectively [ pT| 



To = 



Ao £0 
Ao Do 



(3.17) 



with the renormalized versions of these coupling constants 
being defined as 



_ Z~Zu _ 

u = — Uo Ad p. 

f= J^foAdf^"' 
Z\Zd 

Z-tZd 

T = To , 

ZxZ~ 



(3.18) 
(3.19) 
(3.20) 



whence the upper critical dimension of the dynamical ver- 
tex / comes out to be the same as that of the static vertex 
2, i.e. dc — 4. Also, in equilibrium, Z~ = Z\^ Z~ = Zjj, 
wherefrom we conclude that T = To = 1. 

The analysis is now carried through by considering all 
one-loop Feynman diagrams which contribute to the dif- 
ferent vertex functions listed above. The rcnormalization 
factors are determined in such a way that the renormalized 
vertex functions depend on the renormalized couplings in 
the same way as the zero-loop vertex functions depend 
on the bare couplings at a normalization point, chosen 
here to be (q, cu) — (0, 0) and r = 1 (to one-loop order, 
To ~ /x^), well outside the IR region. The introduction 
of this normalization point (NP) renders the renormal- 
ized couplings dependent on the momentum scale /i, as 
pointed out above. 

From a,2r~^(q,w)|NP and 5q4r~^(q, w)1np, one ob- 
tains the following results for the rcnormalization factors 
Zx and Z^ |22|, 



ZxZr^l 



Zx^l 



I 2(rf-l) 7 Ad 

Uo Ad _e 
d—2-^' ' 
I 2(rf-l) ~ A, 



(3.21) 
(3.22) 



where we have used the value of roc obtained from the 
condition of divergence of the static susceptibility (see be- 
low) and where e — A — d. From these equations one can 
obtain the value of Zt, namely 



Zr = 1 



Uo Ad 



d-2 



(3.23) 



Applying the normalization condition to the derivative 
dq2 r ~~{q^, uj)\-^p , one obtains the following expression for 

Zj/Z^ 

S^l + ^^^M-. (3.24, 

Z^ [a — 2)d lo e 

The normalization conditions for the derivatives of the 
vertex functions 9q2ljj(q, ti;)|Np and 9g2/~j'(q, [x))|np in 
turn yield the following expressions for Zjj and Z~, 



,/o Ad 



A{d + 2) e 



— l^ 



^5"^+4(d + 2) 



1 J rr Ad _ 
Jo To — M 
e 



(3.25) 
(3.26) 
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which become equal in equilibrium, Tq = 1, as they should. 
From the normalization condition for the three point ver- 
tex 5g»r~^j„(-q- p, -w;q,tj;p,0)|NP, one simply ob- 
tains, to one-loop order, that Zg ^ I, which confirms the 
result, valid to all orders, obtained from the condition of 
Galilean invariance. 

On the other hand, from the normalization condition 
for 5(g=p2+p„,2)/~»^^(-q - p, -uj; q, uj; p, 0) \np, one ob- 
tains the following expression for 



Note that in equilibrium this expression is identical to 1, 
i.e., there are no correctio ns to 77 = to one-loop order. 
Eq. ( |3.27D can be used in ( |3.24D to find Z~, 



2(d-l) /o Ad 
{d- 2)d To e ^ 

+ (d-2)d(d + 2)^° tJ e ^ 



(3.28) 



which again reduces to (3.22) in equilibrium. 

Finally, the normalization condition for the relaxation 
vertex a^.r~^^^(-q, -u;; §, |; §, |; §, |)|np yields the 
following expression for the product Z\ Z„ Z^, 



Zx Zi,Z,p = l-'^UQ^fi " 

d-1 J Ad _ 

H ^ Jo — fJ- 

d e 



(3.29) 



from which one obtains, using ( |3^ ) and (|^), the fol- 
lowing result for 



, 3^ Ad 
= I ~ ^uo — fJ- 

4 7 f.l \ ^ - 

(d-2)d(d + 2) ^° \ To) e ^ 



(3.30) 



which is independent of the dynamic coupling /o in equi- 
librium, as it should be. 

The above results for the Z factors are sufficient to pro- 
vide the beta functions which determine the fixed points of 
the theory. However, one still needs to consider the inde- 
pendent renormalizations needed for the susceptibilities, 
since one does not have the usual constraints imposed by 
detailed balance. To define the "static" limit of the intrin- 
sically dynamic model under consideration, we compute 
the response functions for the order parameter and for the 
transverse current, and then take the limit w ^ there. 
In order to obtain these functions, one simply ne eds t o 
take the derivative of the functional J^c, equation (2.17), 
with respect to the external sources /i(x, t) and A(x, t) 
and then set these sources equal to zero. Some formal ma- 
nipulations using the properties of the vertex functions 
yield 

Xo(q,^) = /^o^^(-q, -^r^ 



Xo(q,c^) = r 



ojj 



,(3.31) 



90 r ~ 



[V-vv] 



(-q,-^) ,(3.32) 



respectively [g^|jl0|,|23| . Note that composite-operator ver- 
tex functions enter these expressions, which in general 
implies that new renormalization constants are required 
to remove the UV singularities of the response functions 
(equivalently, one may utilize the Z factors obtained from 
the multiplicative renormalization of the vertex functions 
plus appropriate additive renormalizations p^JIo[ |). Yet 
one may show to all orders in perturbation theory [ p4| 
that 

^ojj + Doq^ - go ^oJ[v^v^] ^) ' (^'^^^ 

and consequently 

Xo(q,c^ = 0) = l , (3.34) 

which means that there is no additional renormalization 
required here. On the other hand, the static limit of the 
order parameter susceptibility is in fact singular, which 
leads us to define the corresponding renormalized response 
function via 

X(q,cj) = Zxo(q,t^) ■ (3.35) 

The new renormalization constant Z is determined by de- 
manding that dq2x~^{<l,uj)\-NP be UV-finite. In this case, 
to one-loop order we obtain the simple result Z — 1, i.e., 
there are no corrections to this order of perturbation the- 
ory. Note that in equilibrium, one has Z^ — Z to all orders 
as a consequence of detailed balance 123] . 

The fluctuation-induced Tc shift is determined from 
the criticality condition Xo^{<l^ 0,uj ^ 0) = at rg = roc 
(tq = 0) with the result 



1 ^ 



1 



kroc 



fc2 



(3.36) 



d-1 Do J 
d Ao 



1 

To 



k^iroc + Do/Xo + k^) ' 



which determines roc implicitly. The momentum integrals 
in (3.36) should be evaluated with a finite upper cutoff, 



which underlines the non- universality of the Tc shift, i.e., 
its dependence on short-distance properties. This expres- 
sion is however sufficient, in its present form, to provide 
for the additive renormalization (mass renormalization) 
necessary to make P~^(q, lu) finite and we have implicitly 



used it to obtain the results (3.21) and (3.22) 



3.2 Discussion of the RG flow equations 

3.2.1 RG equations for the vertex and response functions 

By means of the above renormalization constants, we can 
now write down the RG (Gell-Mann-Low) flow equations 
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for the vertex functions and the dynamic susceptibilities. 
The latter connect the asymptotic theory, where the IR 
singularities become manifest, with a region in parame- 
ter space where the loop integrals are finite and ordinary 
'naive' perturbation expansion is applicable, and follow 
from the simple observation that the 'bare' vertex func- 
tions do not depend on the renormalization scale yit, 



and 



(3.37) 



with {flo} — Aq, Ao, J^o, Dq, To, uq and g p . Replacing the 
bare parameters and fields in Eq. (3.37) with the renor- 
malized ones, we thus find the following partial differential 
equations for the renormalized vertex functions. 



{a} 

xr^~„., (Ai, {q,w};{a}) = . 



Here, we have introduced Wilson's flow functions 

d 



In Z, 



4, 



h-iZ 



and 



In — 

"^0 



(3.38) 

(3.39) 
(3.40) 

(3.41) 



(the index '0' indicates that the renormalized fields and 
parameters are to be expressed in terms of their bare 
counterparts prior to performing the derivatives with re- 
spect to the momentum scale fi). Note that (~ = (j =0, 

(g = d/2 — 3 as a consequence of E qs. (3.6 ) and (|3.14D . 

The Gell-Mann-Low equation ( 3.38 ) is readily solved 
with the method of characteristics /i —^ fi£] this defines 
running couplings as the solutions to the first-order differ- 
ential RG flow equations 



da{£) 



Ui)a{e), a(l) =a 



(3.42) 



The solutions of the partial differential equations ( 3.38 ) 
then read 



(3.43) 



exp < t: 



xr^~j,^jAi^,{q,c.};{a(f)}) . 

In the same manner, one can solve the RG equations for 
the dynamic susceptibilities, with the results 



^(m, {q,^};{a}) = 

= X{tii,{(i,uj};{a{£)}) 



(3.44) 



x(m, {q,^};{a}) = 

= exp < 



(3.45) 



where, in analogy with Eq. ( 3.3£ ), 

InZ 



(3.46) 



In terms of th e ren or mali zed couplings u, / and T, as 
defined by Eqs. ( 3.18 )-( ^.2C| ), we find, using the results 
of subsection 3J for the Z factors, the following explicit 
results for the zeta functions to one- loop order: 



■l-T 



{d-2)d{d + 2) 

i / 

V {d - 2)d{d + 2) 
C- + 0(^2,?), 



T ' 
l-T 
T 



Cd 

Cr = 



2{d-l) 



f, 



{d-2)d 

2(rf-l) / 
{d-2)d T 

J 

4(d-f2) ' 
1 



{d~2)d{d + 2) 



l-T 
T 



4(d + 2) 



d-2 
3 ^ 



— ~f -\- — ?/ — 

2 {d-2)d{d 



l-T 



T 



(3.47) 

(3.48) 

(3.49) 
(3.50) 

,(3.51) 

(3.52) 
(3.53) 
(3.54) 
(3.55) 



with Cg = d/2 — 3, identically. These results enable us 
now to study the scaling behavior of the non-equilibrium 
model H with dynamical noise in the vicinity of the differ- 
ent RG fixed points, which are given by the zeros of the 
appropriate RG beta functions ({«} = u,f and T) 



d_ 
9/i 



According to 



de 



(3.56) 



(3.57) 



these govern the flow of the effective couplings u, /, and T 
under scale transformations n — > and the fixed points 
{v*} where all f3v{{v*}) = thus describe scale-invariant 
regimes. 

The RG analysis of the theory then requires the study 
of the behavior of three independent coupling con stant s 
under the RG flow, the static coupling u defined in ( |3.18 ), 
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the mode-coupling vert ex / defined in (3.19), and the pa- 
rameter T defined in ( 3.2C ), which denotes the ratio of 
the noise temperature of the order parameter to the noise 
temperature of the transverse mass current, with the cor- 
responding beta functions being f3~, /3~and Pt- 

Evaluating the Gell-Mann-Low flow equations near an 
IR-stable RG fixed point, we may derive the following scal- 
ing laws for the two-point correlation functions of the or- 
der parameter and conserved currents, respectively, 

C^(r,q,^)=,-2--^+''^'(7,,(^-^,^) , (3.58) 



and for the order parameter susceptibility 

T LO 



ql/u ' 



(3.59) 



(3.60) 



Here, the different critical exponents are given in terms of 
the following fixed-point values (indicated by a '*') of the 
zeta functions 



2{d^ + d-A) 



[d- 



2)d{d^ 
T 

"4(d+2) 



2) 



(3.70) 



where we have used Eqs. (3.50)^( 3.55 ) for the one-loop 
Wilson zeta functions. From Eq. ( 3.701) one sees that there 
exists a fixed point either if T = 1, which is the ordinary 
equilibrium fixed point at which the noise temperatures 
become equal, or when the mode coupling vertex / = 0, 
which can yield a non-equilibrium fixed point (see below). 

As a prelude to the study of the behavior of the model 
at the non-equilibrium fixed point and in order to render 
the analysis more transparent, we start by reviewing the 
results for the eq uilibri um fixed point. Firstly, we notice 
that in equation ( 3.70 ) one has /3t > for T > 1 and 
/3r < for T < 1. This implies that in the IR regime 
£ — > 0, r increases if T < 1 and it decreases if T > 1. 
The equilibrium fixed point is thus stable with respect 
to perturbations in the noise temperature. Secondly, at 
this fixed point one obtains the following values for the 
coupling constants u* and /*, 





(3.61) 




(3.62) 


= ~C , 


(3.63) 


= 4 + Ca , 


(3.64) 


^V = 4 + C~, 


(3.65) 




(3.66) 


^J=2 + C5. 


(3.67) 



The first three exponents correspond in equilibrium to the 
static critical exponents (with 77 = 77^ ) , while the last four 
yield the dynamical critical exponents (again, in equilib- 
rium = and zj = zj ) . 



3.2.2 RG fixed points and critical exponents 

The relevant RG beta functions, namely /?-, /?r and /3t, 
are given, to one-loop order, by the following expressions: 



-uU-lu+'^^-^f], (3.68) 



2 {d-2)d T 



and 



/37=-/(e + 2CA + Cz?-Cx 

_ -/ nd^ + Ud - 48 ~ 
^ \ 4:{d ~ 2)d{d + 2) ■' 

2{d^ + d-A) P 
{d-2)d{d+2) T 



/3t = T(Cx + Cd-Ca-C5 



(3.69) 



19 



O(e^), 

f 0(6^) 



(3.71) 
(3.72) 



where we h ave p erform ed a n expansion around d = 4 in 
expressions ( 3.68 ) and ( 3.69 ) and where we have used the 
fact that (3^ is independent of / when T = 1. When these 
values are s ubsti tuted in the expressions for the zeta func- 
tions ( |l47| ) to (|l55| ) , one obtains for the critical expo- 
nents, Eqs. ( |3.61 ) to ( [3.67| ), the results 



18 



^^=4-— e + 0(e'). 



e 

19 



+ 0{e^). 



(3.73) 
(3.74) 

(3.75) 
(3.76) 



The values for the static critical exponents are (due to 
the existence of detailed balance) the ones one obtains if 
one performs an RG _analysis of the static model described 
by the free energy (2.1). As for the dynamical exponents, 
whereas zj shows a small negative correction due to the dy- 
namical coupling between the two modes, the dynamical 
exponent z^, which describes the behavior of the charac- 
teristic frequency of the order parameter (q) oc q^^ , dis- 
plays a strong suppression already in d = 3 (e = 1). This 
can be understood from the fact that in the absence of dy- 
namical coupling and at the critical point, the dynamics 
of the conserved current is purely diffusive and therefore 
much faster than the dynamics of the also conserved order 
parameter which shows the characteristic critical slowing 
down [ p6[ . The existence of dynamical coupling between 
the two modes will then provide the slower order param- 
eter field with an additional and faster decay channel, so 
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that < 4. Notice furthermore that Eq. (3.69) imphes for 
any fixed point with non-trivial and finite mode-coupling 
< /* < CXD, in equilibrium (where Ca = C'j) the relation 
C\ + Cd = — e = — 4, i.e., the exponent identity 



z^ + =6 



d+2 



(3.77) 



which is of cours e satisfied by the explicit one-loop results 
( ^J5|) and ( |376|) . 

If one now takes / = and T finite in Eqs. ( 3.68| ) to 

( 3.70 ), one obtains a line of fixed points (i.e., one fixed 
point for each value of T), provided that u = | e -I- O(e^). 
This line of fixed points corresponds to a system where 
the two modes are completely decoupled, as can be seen 
from the computation of the critical exponents. One has 



z^ = z^ = A + O(e^) , 



2J 



(3.78) 

(3.79) 

(3.80) 
(3.81) 



with the last equality holding identically, given that / = 0. 
The static critical exponents and the dynamical exponents 
z^, z^ are simply the exponents one would obtain for the 
equilibrium purely relaxational critical dynamics model B 
with conserved order parameter, with the O(e^) term rep- 
resenting two-loop contributions coming from the static 
vertex. For the dynamic critical exponent, one knows that 
in fact 

Z.0 = 4 - (3.82) 



exactly |l|,[|. However, from Eq. ( p9| ) we see that the 
model B fixed point is IR-unstable against perturbations 
in / since Pj< for small /. 

The cases f — with T = or T = 00 require extra 
care, as one needs to consider the flow of /' = //T in the 
first case, and the fiow of / = / T in the second case, since 
these two quantities might, respectively, have a finite value 
at the hypothetical fixed points. If / = and T = 0, then 
/3p = — e/' with /' = f /T, which shows that a fixed point 
cannot exist for any finite value of /'. Physically, the ab- 
sence of a fixed point is again due to the fact that the order 
parameter displays critical slowing down compared to its 
diffusive relaxation away from the critical point. The con- 
served current is, on the other hand, always governed by 
faster, diffusive dynamics. If the conserved current were to 
be placed at infinite temperature, it would be slaving the 
order parameter, which is impossible. Thus no fixed point 
corresponding to such a situation exists, in contrast with 
the isotropic non-equilibrium SSS model (lO|. There, the 
order parameter is not conserved, thus permitting a T = 
fixed point (albeit an unstable one) where the purely dif- 
fusive (as there is no feedback from the order parameter) 
conserved variable actually becomes the slower dynamical 
mode, with z^ = 2 — e/2 + O(e^) for the order parameter 
011 ■ 



If / = and T = 00, one finds, on the other hand, for 
the RG beta function corresponding to the new coupling 

f = fT: 

^/--G-4(dT^^")^"' ^'-''^ 

which has a non- equilibrium fixed point for /* = 4((i + 2)e. 
It is also easy to see that, since /3jr<Ofor/</*,/3j>0 
for / > /*, this genuine non-equilibrium fixed point is 
IK- stable with respect to perturbations on the value of /. 
However, it is unstable against perturbations on the value 
of T, since /3t > for large T. The values of the critical 
exponents at this unstable fixed point are given by 



Z4, 



+ 0(6^) 

O(e^) , 



z, = 2 - e + 0(6^) = d~2 + 0{e^) . 



(3.84) 

(3.85) 

(3.86) 
(3.87) 
(3.88) 



Eq. ( |3.87 ) holds to all orders given the structure of the 
vertices at this fixed point. The O(e^) corrections to the 
static critical exponents and to z^ — z^ stem from the 
two-loop contributions induced by the s tatic vertex only. 
We thus expect the exact relation ( 3.82| ) to hold here as 
well. Moreover, as we have C.x — C.~ at T = 00, the very 

existence of any non-trivial, finite fixed point < / < 00 
demands that -|- ^~ = — e, i.e., z^ + Zj = 6 — e~d + 2, 
which ge nera lizes the corresponding equilibrium scaling 
relation ( 3.77 ). With C\ = —77^ we therefore arrive at 



z^ ^ d - 2 + rj^, 



(3.89) 



One sees that this fixed point is characterized by model 
B exponents for the order parameter and by anomalous 
noise correlations for the conserved current, giving rise to 
a value of Zj ^2 for d < dc — 4. Again, the model B expo- 
nents obtained for the order parameter can be understood 
from the fact that the conserved current is effectively at 
zero temperature, and therefore its dynamics does not in- 
fiuence the dynamics of the order parameter, being how- 
ever affected by it through the residual one-way coupling 
between these two dynamic variables. Such a fixedpoint is 
also present in the non-equilibrium SSS model with 
similar model A-type behavior for the order parameter, 
and anomalous noise properties for the dynamically cou- 
pled conserved fields. In the present context its existence 
is due to the fact that here the slower variable is at higher 
temperature, slaving the faster conserved modes, and suf- 
fering no feedback from the latter. The characteristics of 
this novel, unstable non-equilibrium fixed point should be 
contrasted with the pure model B fixed point discussed 
above where there is no coupling at all: f ~ 0, and where 
the temperature ratio T, albeit finite, does not matter, for 
the two Langevin equations are fully decoupled. The RG 
fixe d po int structure, and their stability is summarized in 
Fig. |t. 
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1/2- 



noneq. model H 

-o 



, eq. model H 









model B 



1/2 



1 



Fig. 3.1. Equilibrium and non-equilibrium fixed points for the 
isotropic non-equilibrium model H, plotted for d = 3 (e = 1). 
The equilibrium fixed point in the center of the flow diagram is 
infrared-stable. The model B fixed line is unstable with respect 
to the mode-coupling, the non-equilibrium fixed point at T = 
oo is unstable in the T direction. 



4 The anisotropic non-equilibrium model H 

In this section, we study the critical behavior of our non- 
equilibrium version fo r m odel H w ith dynamical noise, as 
defined through E qs. ( 2.4) and (2.6) and by the anisotropic 
noise correlators ( 2.11 ) and ( 2.12 ). We start by computing 
the Tc shift from the static susceptibility. As a consequence 
of the spatially anisotropic conserved noise with Tq < Tg , 
it turns out that the transverse momentum space sector 
with lower noise temperature softens first. Thus, at the 
critical point, the longitudinal sector remains uncritical 
('stiff'), similar to equilibrium anisotropic elastic p hase 
I and at Lifshitz points [ 



4.2 



transitions |^ and at Lifshitz points In Sec 
we turn to the perturbational rcnormalization of the two- 
temperature non-equilibrium model J to one-loop order 
and finally discuss the resulting RG flow equations. 



4.1 Dynamic field theory and the anisotropic Tc shift 

The dynamic field theory which describes the anisotropic 
model H has already been present ed in Sec. ^, th e dynamic 
functio nal J bein g give n by Eqs. ( 2.14 ) to ( 2.16| ), with the 
choice ( ^.llD and ( |2.12| ) for the noise correlators. However, 
the computation of the shift in the critical temperature 
shows that, as it stands, this theory is not fully consistent. 
The shift in the critical temperature is determined, as in 
the isotropic theory, by the condition Xo^i^i^) = in 
the limit q — > 0, w — > 0, where xo(q, w) is the order 
parameter dynamic susceptibility, whic h is gi ven in terms 
of the relevant vertex functions by Eq. ( 3.31 ). 

From the ensuing expression (to one-loop order), we 
may determine the fluctuation-induced shift of the criti- 
cal temperature. Because of the dynam ic anisotropy ap- 
pearing in the noise correlators (2.11) and (2.12), how- 
ever, the result depends on how the limit q — > is taken; 
upon defining g|| = qcosO and q±^ = 5 sin 6* and with 



(A^/Ao)(I?o/i^o^) denot- 
ing the temperature ratios in the different sectors, we find 



roc{0) -~ 



d Do 
1) 



1 



1 



XUd + 2) 

d Do 

9lid+l) 
XUd+2) 



(tS - 1) 



k roc + 
1 



(4.1) 



k^{roc + Do/Xo + k^) 
1 



k roc + fc^ 



(To^ 



1 



fe k^{roc + Do/Xo + k^) 
1 



Xld{d + 2) k^roc + Do/Xo + k^) 




cos^ + 










\Xo 


Do J 





in contrast with Eq. ( |3.36 ) for the isotropic model H. Here, 
(i|l and d± are, respectively, the dimensions of the parallel 
and transverse subspaces, with d= d\\+d±. As Tc = T° -|- 
roc, the phase transition will occur at the maximum of the 



{ Ao^ 




< 






I Ao 


Do ) 


i Ao 


Do j 



or. 



equivalently, Tq^ - 1 < {T^ - 1)^, occurs at 61 = n/2. 

^0 

The d^-dimensional transverse sector in momentum space 
thus softens first, and the true Tc shift is given by 



V2/ d Do \ 2 Jk^Oi 



1 



1 



fc2 



(4.2) 



alid+l) \\ 
'Xl{d + 2r^ ^ 

di D. 



1 



d Do 
X I (i+ 1 



2 .Ik roc 



k k^{roc + Do/Xo + k^) 
1 



2 
.95 



2 

dl 



1 



A^(d + 2) 



(To^- 



Ik k^{roc + Do/Xo + k^)^ 

where again this non- universal quantity must be computed 
with a finite UV cutoff. For T^ = T^ = Tq, we recover the 



result ( 3.36 ) for the isotropic model. Notice that dynam- 
ical anisotropy (Tq =/= TJ^), combined with the reversible 
mode-coupling terms, has a very drastic effect here: It ren- 
ders the system soft only in the momentum subspace with 
lower effective noise temperature. This effect has a sim- 
ple physical interpretation: The Tc shift is due to thermal 
fluctuations, which are reduced in the transverse sector 
(T(^ < Tq), and therefore lead to a comparatively stronger 
downwards shift in the longitudinal sector. This result is 
completely analogous to our earlier flndings for the non- 
equilibrium model J (describing isotropic ferromagnets) 
with anisotropic noise [ll|| . 

In order to characterize the critical properties of our 
model, we may neglect terms oc gjj in the stiff momen- 
tum space sector, because = — roc{0 = 0) remains 
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positive at the phase transition where t^^- ~ tq — rQc{0 — 
7r/2) vanishes. In analogy with the situation at anisotropic 
elastic structural phase transitions ||2^ , or with Lifshitz 
points in magnetic systems with competing interactions 
pSf , as well as driven diffusive systems ||^, we thus have 
to scale the soft and stiff wave vector components differ- 



/i, whereas [q\\] 

,0 



Consequently, 

4 



ently, [q± 

while [Xq] = fi^ and [Dq] = h'^, if we choose [uj\ = we 
find for the longitudinal noise strengths the scaling dimen- 
sions [Xq] = and [ZJq] = which implies that they 
become irrelevant under scale transformations. Allowing 
for distinct couplings in the different sectors, one finds in 
the same manner that the ratios [Aq/AJ-] = [Dq/Z?,^] = 
[Ajjuil/A^u^] = aj^d [gl/go] = aU have nega- 
tive scaling dimension. Thus, for an investigation of the 
asymptotic critical behavior, the longitudinal parameters 
may be neglected as compared to their transverse counter- 
parts, and can all be set to zero in the effective dynamic 
functional. ^ 

Upon rescaling the fields according to V-'o ^ [Xq /Xq)^/'^ 

^0, V^o^- (A(f/A^) 1/2^^0, }d - (D^/Dono and ^ 
(D^/D^)^/^ and defining 



Co = 



go = 




Uo = 



^0 



and omitting the labels '_L' again for Aq and tq, the ensuing 
effective Langevin equations of motion become 



dtpo 

at 



= An 



coVi|+Vi(ro-Vi; 



^0 



(4.4) 



and 



^=r [i^oVijo + go'^^Mro - Vi)^o + C] , (4.5) 

where the transverse projector in the soft (±) subspace 
T is given, in Fourier space, by the expression T"''(k) = 
_ k1k^/k'j_. The noise correlators in turn read 



(r/(k, i/)?7(k', t.')) - 2Ao kl S{k + k') S{i^ + v') (4.6) 



and 



(C"(k, i.)C''(k', v')) = 2D^ kl <5(k + k') 5{v 



(4.7) 



where again, for convenience, we have used the Fourier 
space representation. These equations define the two-tem- 
perature non- equilibrium model H. In order to perform the 
RG analysis, one represents these equations in the form of 
a dynamic functional, precisely as in Sec. H above. 



We emph asiz e the fact that the anisotropy of the Tc 
shift in Eq. (4.1) only occurs in the contribution oc g^, 
i.e., the anisotropy in the Tc shift is due to the purely 
dynamical mode-coupling terms. In the non-equilibrium 
model B with dynamical anisotropy , the criticality con- 
dition for the response function remains isotropic, at least 
to one-loop order. Thus, if one does not assume differ- 
ent critical temperatures in the purely diffusive non-linear 
Langevin equation to begin with, these are not generated, 
and one is not immediately led to the two-temperature 
model B as the correct effective theory for the phase tran- 
sition. In the presence of reversible mode-couplings, how- 
ever, anisotropic noise correlations, specifically for a con- 
served order parameter field, have a much more drastic 
effect: For both models J [|ll| and model H such viola- 
tions of the detailed-balance conditions render the system 
inherently anisotropic at criticality, and certainly prevent 
any restoration of the equilibrium critical properties. 



4.2 Renormalization of the two-temperature model H 

We start by noticing that, as with the isotropic non-equi- 
librium model H, the two-temperature model H, being a 
genuinely non-equilibrium model as well, does not allow us 
to invoke a fiuctuation-dissipation theorem in order to re- 
late vertex and response function renormalizations, and we 
have to compute almost all the Z factors independently. 
These consist of the wave function renormalization fac- 
tors Z^, 



and the 



"^J' ^ ^'^'^^'^ ^ 
coupling constant renormalizations, which we define here 

through 



A 


= ZxXo , 




(4.8) 


c 


= ZcCa, 




(4.9) 


D 


= Zd Dq , 




(4.10) 


T 


= Zt- To /i^^ with 


To = ro- roc , 


(4.11) 


U 


= Z~uoA{d\\,di^) 




(4.12) 


g 


= Zy^g„A{d\\,d, 


)l/2^(d+d||)/2-3 ^ 


(4.13) 


g 


= Z^'^ go A{d\\ , dj 


)l/2^(d+d||)/2-3 ^ 


(4.14) 



where A{d\\, d±) is given by 

r(3-d/2-d||/2)r(d/2) 



A(d||,rfj 



2d-iTTd/^r{d±/2) 



Again, this factor is included because it appears 
quent formulas. 



(4.15) 
in subse- 



As before, Eq. (3.1) implies that the relations 
(3.5) for the wave function renormalizations hold, 
more, the vertex structure of the model leads to 



(|3J) and 
Further- 



r,^,^ (q||, - 0, ^ = 0) = Ao CO q^ , (4.16) 
which must hold to all orders. This entails that 



ZaZ, = 1, (4.17) 
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which leaves us with one Z factor less to determine. One 
has thus a total of seven independent renormalization fac- 
tors. As the vertex functions or their derivatives which 
we must render finite are the same as for the isotropic 
model H, given at the beginning of Sec. 



3.1, with the ex- 



ception that one has to substitute the derivatives dq with 
respect to q by derivatives dq^ with respect to q^, one 
sees that the renormalization conditions on these vertex 
functions determine all the renormalization factors, i.e., 
there is no freedom left, as opposed to the isotropic model 
H, to arbitrarily fix one of the renormalizations. To these 
renormalization factors, one ad ds, as above, the renormal- 
ization factor Z, defined in Eq. ( 3.35| ), necessary to render 
the order parameter susceptibility finite. It is determined 
by the condition that 9^2^x~^(fl7 '^)|np be UV-finite |^ . 
Furthermore, Galilean invariance still holds in the two- 
temperature model (see Appendix which entails that 

(4.18) 



to all orders, although, as pointed above, one can now no 
longer suppose that — 1. 

Next we introduce the coupling constants 



fo = 



Up 

goffo 
\ n /II /2 ' 



\ n /II 



and their renormalized counterparts 

z'i||/2 



V = Z~Zs 

u A 



/ 



,l-d||/2 



■/o A(d||,d_L)^~ 



Zd 



/ = 



{z,z~Y'' 



foA{dpdj_)n- 



Z~ 

Z. Zd 



(4.19) 
(4.20) 
(4.21) 
(4.22) 

(4.23) 
(4.24) 

(4.25) 
(4.26) 



where e = 4 — c? — . Notice that all these renormalized 
non-linear couplings become scale-invariant at the reduced 
upper critical dimension 



rfc(rf||) =4-d|| 



(4.27) 



Such a lowering of the critical dimension is typical of mod- 
els with anisotropic scaling [ p7[|2^ ,|5t, for the fluctuations 
are critical merely in the transverse sector. Employing the 
renormalization conditions, one then obtains to one- loop 
order the following results for the Z factors 



6(d^-2)(rf^-l) - A(d||,d^) _ 
^'-^ + Jd^2jdKdZT2)^" -e ^ 



(4.28) 



Zr = l 



Zd = 1 



vq dj_-l - \ A{dii,dj_) _ 
fo M 



d + dn -2 V 2 



6(di-2)(rfx-l) J Md\\,d±) _ 
Jo M 



{d-2)d^{dj_ + 2) 
1 



4(rfi + 2) 



A(d||,dj.) _ 
70 M 



2{d^ - 1) A(d||,d^) 

id~2)d^^' -e ^ 

6(d^-2)(rfx-l) A{d^^,d^) _ 
(d - 2)d^{d^ + 2) ■'° e ^ 
1 , A(dn,d_L) 



(4.29) 
(4.30) 
(4.31) 



^J-^+4(rf^ + 2) 



fo- 



4d_L(d_L + 2) 



~ A{d\\,di_) _ 
fo M 



3 A(d||,d^) _^ 
Z~=l--vo M 

2{d^ - 1) A(d||,dx) 

+ {d-2)d^-^' e ^ 

d^-1 / 12(di-2) \~ Aid\\,d±) _ 

i - — TT I Jo M 



(4.32) 
(4.33) 



(d-2)(di +2) 



Z~ = 1^ 

a {d-2)di_ 

4 



A(d||,d^) _ 
/o A* 



d_L(d± + 2) 

J A{d\\,di_) _ 
X /o M 



3(di - 1) 



d^+2j e 
d^-2 rf^ + 16 



d-2 



16 



4d_L(d± + 2; 



~ A{d\\,d^) 
■ Jo A* 



(4.34) 



with Zg — Z-^ ^ , a result which can be confirmed explicitly 
to one-loop order from the renormalization condition of 
dqiF-^^^l-q - p,-u;;q,uj;p,0) \np- Subsequently, ren- 
dering 9^2^x~^(q, — 0)|np UV-finite yields the addi- 
tional Z factor for the response function 

4(dx-4)(d^-l) A(d||,dJ 
^-^+(d-2)dx(d^ + 2)^" ^ ^ -^^-^^^ 

These Z factors can now be used to compute the relevant 
beta functions and Wilson zeta functions of the theory. 



4.3 Discussion of the RG flow equations 

Being in possession of the expressions for the Z factors to 
one-loop order, one can, i n an analogous manner to what 
was done in sections ^.2.1 and 3.2.2 , compute the relevant 
Wilson zeta functions, also to one-loop order. These are 
given by 



Cj = - 



{d - 2)di_ 
1 



6 



4(di 



(rf±-2)(dx-l) 
[d-2)d±{d±+2.) 



Ad^{di 



(4.36) 
(4.37) 
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c = 

Cd 



4(dx~4)(d^-l) 
'(d-2)di(di + 2) 
6(d^-2)(rf^-l) 
(d - 2)dj_{d± + 2) 



f, 

f: 



Cr = -2 



4(d_L + 2) 
1 



2(rf± - 1) 



d + d\\-2 
1 3(dj 



d + d||-2 {d~2){d± 
_ 3 2(rf^ - 1) 
2 (d - 2)di 



(4.38) 
(4.39) 
(4.40) 
(4.41) 

(4.42) 



rfi - 1 



1 - 



2C~ = " "II 



d + dn-6 
4 



12(d^ - 2) 

{d-2)(d^+2) 
4 



c!x(d± 2) 
d 

Ad^{d^+2) 



{d - 2)dj_ 
3{dj_-2){dj_ - 1 



d, - 1 



d-2 



d_L + 2 
(i_L + 16 
16 



/ 
/ 

(4.43) 



with Cc = -Ca from Eq. ( |4.17D and 2Cg = d + d|| - 6 - Cj 
from (4.18). These equations reduce to the equihbrium 
ones in the limit d\\ = 0, i.e., d± = d. 

In the anis ot ropic two-temperature model H, the scal- 
ing laws ( 3.58| )-(3.60) generalize to 



Cj(T,q||,q±,tj) 



911 



,9 1 1± 1± 



-2-zi 



= 1l 



-2+J7 - 
9^ X 



II V ' l + ' zj 



and the critical exponents are defined via 



(4.44) 
(4.45) 

(4.46) 





(4.47) 




(4.48) 


= -c , 


(4.49) 


Zi = l-^ = l + ^, 
2 2 ' 


(4.50) 


= 4 + a , 


(4.51) 


= 2 + , 


(4.52) 



where the exponent A originates from the intrinsic aniso- 
tropy of the system |^,^,^ . 



From the zeta functions ( 1.3£ )-( 4.43 ), one can com- 
pute the beta functions for the coupling constants f and 
/, which determine the fixed points, with / being given 



hy f — \/ ff. To one-loop order, these beta functions read 



3 2{d^ - 1) 
2 " (d - 2)d_L 



(4.53) 



di 



d± 



1 ^ 3(d|| -4)(rf^-2) ^ ^ 



(d-2)(d_L + 2) 



2(Cg + 1) - + 



1 Ca 



f = -ef 



1 



2{di_ + 2) 
d 

'AdJdZ+Y) 



3(d|| - 2){di_ - 2)(di - 1) 
{d-2)d^{djL + 2) 



^7 = 



2(C~+1)-Cd+ -1 



dl - 1 - 



/ 

(4.54) 
(4.55) 



(d - 2)dA 



f 



1 



d^{d^ + 2) 
d 



3(d|| -6)(c!i-2)(di-l) 



d-2 



4di(di + 2) 
and it follows from / = \fjf that 



/3f 



/ / 



(4.56) 



which is used, together with (4.54) and (4.55), to compute 

Since the beta functions (3f and /3~ do not depend on 
the static coupling u, the determination of the fixed points 
reduces to the solution of the system of quadratic equa- 
tions given by /?/ — j3j— 0, the equation for I3j being au- 
tomatically satisfied. Introducing the following condensed 
notation 



2 (d - 2) d^ {di_ + 2) 



7 = 



{d - 2)d± - 6(d|| - 2){dj_ - 2){di_ - 1) ' 

4(d-2)d_L {djL + 2) 
(d- 2) (di - 8) - 12 (d|| - 4) (di - 2) (di - 1) 
{d-2)d 

2[{d - 2)di_ - 6(d|| - 2){di_ - 2){di_ - 1)] 

8[(dj. -l)(d^ + 2)-2] 
(d-2)(di-8)-12(d||-4)(di-2)(di-l) 



(4.57) 
, (4.58) 
, (4.59) 
, (4.60) 



one can write the solutions of the system of quadratic 
equations as 



a- a- 2a7c ± [a - a)^ + Aacaj 

f± = TTf T\ 4-61) 

27 7c- 1) 



~ _ a - a - 2ac7 ± ^/ {a - + \aca^ 

^± - 2c(7c-l) ^^-^^^ 
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The existence of fixed points depends on the existence of 
at least one positive root for each of these two equations, 
which gives a series of conditions on the coefficients a, c , 
a an d 7. H owever, it is immediate to see from Eqs. ( 4.61 ) 
and ( |4.62| ) that, even when those solutions exist, they di- 
verge if 7c = 1. Taking into account that d = d_i_ + d\\, this 
condition defines a sixth-order equation determining the 
relation between d± and . With the minimal subtraction 
prescription that d = A — d\\, i.e., setting e = in this equa- 
tion, one obtains the numerical solution d|| = 0.838454; 
already for < dy =1, the RG flow takes the mode cou- 
pling to infinity. This result is completely analogous to 
the result we have previously obtained in our study of the 
two-temperature non-equilibrium model J p| . 

Formally, and following our study of model J, we may 
expand about the equilibrium model H, and thus obtain 
critical exponents in the limit <C 1. To first order in 
due, we find 



« 24 
f = — e 
•' 19 

7* 24 

^ ^19^ 



1442 , 

11246 

I 1 

6137 

143 
323 



d||e. 



leading to the critical exponents 



21 

'^ = -19^11^' '^'^^-646 ^11" 

-1^1 41 , 

ly =2 e diie , 

3 646 " ' 

, 18 2820 , 

z = 4 e due , 

19 6137 " ' 

„ 1 372 

zi = 2 e - 



19 



6137 



dye. 



A 1 9 1410 , 

A = l e due. 

19 6137 " 



(4.63) 
(4.64) 
(4.65) 

(4.66) 
(4.67) 
(4.68) 
(4.69) 
(4.70) 



Notice that this procedure amounts to an expansion 
with respect to two dimensional parameters, namely e = 
4 — d — d|| , and dye. Moreover, the divergence of the non- 
expanded fixed point /* at d|| « 0.838 454 i ndic ates t hat 
an extrapolation of the formal results ( |4.66 ) to ( 4.70| ) to 
any physical dimension d|| > 1 is unlikely to work. On 
the other hand, we cannot exclude that, also for model 
H, this divergence merely represents an one-loop artifact, 
and is cured if one calculates the RG beta functions to 
higher loop orders. Yet another possibility might well be 
that the divergence of /*, /* and v* indicates the absence 
of a simple non-equilibrium stationary state of the two- 
temperature model H in the vicinity of its critical point. 
For example, in a uniformly rather than randomly driven 
non-equilibrium version of model J, a similar divergence 
has been found recently |]30| . In that case, computer sim- 
ulations have revealed that the system enters a regime 
of spatio-temporal chaos at long times; perhaps the ab- 
sence of a finite RG fixed point in the randomly driven 
two-tcmperature models J and H might indicate similar 
behavior. A somewhat less drastic implication may be 



that merely perturbation theory breaks down, and non- 
perturbative approaches could possibly characterize the 
scaling behavior at the transition of the two-temperature 
model H successfully. 



5 Summary and final remarks 

We have studied two non-equilibrium generalizations of 
the dynamical model H with both conserved scalar order 
parameter and dynamically coupled conserved transverse 
currents that describes second-order liquid-gas or binary- 
fluid phase transitions. Specifically, we were interested in 
the effect of detailed balance violations on the asymptotic 
critical behavior. We have investigated both (a) isotropic 
violations of the equilibrium conditions, which can be for- 
mulated in terms of different effective noise temperatures 
for the order parameter and conserved currents, respec- 
tively, and (b) spatially anisotropic detailed balance vio- 
lations, i.e., dynamical noise which is governed by different 
strenghts in longitudinal and transverse momentum space 
sectors. 

In principle, there are several possible scenarios: (1) 
In the vicinity of the critical point, detailed balance may 
effectively become restored as a consequence of the di- 
verging correlation length that essentially averages over 
the different local noise temperatures; (2) a novel, stable 
renormalization group fixed point may emerge that de- 
scribes a new universality class with genuine non-equilibri- 
um scaling behavior; (3) there might be no stable RG fixed 
point at all, indicating perhaps complex spatio-temporal 
chaotic behavior rather than a simple stationary non-equi- 
librium state. 

We find that situation (1) applies to the isotropic non- 
equilibrium model H, while scenario (3) appears to de- 
scribe the effective two-temperature model that emerges 
upon allowing for spatially anisotropic noise correlations. 
Remarkably, case (2) is never realized in any dynamical 
model with reversible mode-couplings. In fact, a surpris- 
ingly simple overall picture emerges (we have already pre- 
sented a brief overview in Ref. ||T2| ) . Namely, quite gener- 
ally, the equilibrium dynamical models as listed in Ref. Q] 
with non-conserved order parameter turn out to be quite 
robust against detailed-balance violations. The purely re- 
laxational models A and C do not even have a genuine 
non-equilibrium fixed point at all. For the SSS model, 
generalizing models E (for planar ferromagnets, n = 2) 
and model G (for isotropic antiferromagnets, n = 3) to 
arbitrary order parameter space dimension n, two non- 
equilibrium fixed points do exist, corresponding to ra- 
tios T = and T — 00 for the noise temperatures of 
the order parameter and dynamically coupled conserved 
fields, but neither of these is stable. Thus, near the criti- 
cal point, detailed balance becomes eventually restored, 
and the asymptotic critical exponents are those of the 
equilibrium model [Q. Of course, such systems might re- 
main in the crossover region for quite a while, masking 
the asymptotic regime. Essentially, this scenario (1) also 
applies when the conserved noise for the coupled non- 
critical fields is rendered anisotropic as well. Additional 
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fixed points emerge, but the isotropic equilibrium one re- 
mams stable [0. 

When detailed balance is violated isotropically in the 
models B, D (purely relaxational) , J (isotropic ferromag- 
nets, n = 3) and H with conserved order parameter, basi- 
cally the same statements apply, and scenario (1) is real- 
ized again [ p^|Jl^ . However, once one allows for spatially 
anisotropic or dynamical noise, separating a soft trans- 
verse and stiff longitudinal momentum space sector, which 
enforces anisotropic scaling, the behavior changes dramat- 
ically. In the relaxational two-temperature models B and 
D, the ensuing asymptotic theory however turns out to 
be equivalent to an equilibrium model with long-range 
correlations of the uniaxial dipolar or elastic type ||9|,[l2|. 
This corresponds to case (2) above. In stark contrast, 
in the anisotropic non-equilibrium versions of models J 
and H which are characterized by relevant reversible 
mode-couplings to additional conserved variables, no sta- 
ble renormalization group fixed point can be found (at 
least to one-loop order), which represents scenario (3). We 
do at this point not really know what the absence of an RG 
fixed point means physically in this situation; perhaps, as 
in the uniformly driven non-equilibrium model J [pG] , the 
long-time behavior is governed by spatio-temporal chaos. 
It would certainly be worthwhile to explore this issue fur- 
ther, e.g., through computer simulations. 
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A Galilean invariance in the dynamic model H 



In this appendix we derive the basic Ward identity which 
was used in the analysis of the non-equilibri um i sotropic 
and two-temperature model H, namely Eq. (3.14), which 
follows from the Galilean invariance of the Langevin equa- 
tions which describe the model [ p^ . In order to prove such 
an identity, we write the source-free equations ( ^.4[ ) and 
(2.6) in a slightly different form: 



and 



dt 



dip 

dt 



+ 50 Jo ■ VV'o = AoV2(ro - 



V'^o+V, (A.l) 



.9oJo • Vjo 



= r [i^oV^jo + goVMra ~ + C] , (A.2) 

where we have added the convective term .go Jo ' ^Jo 



eluded explicitly in the analysis, as it generates diagram- 
matic contributions proportional to g^/ Dq. Since the scal- 
ing dimension of this effective coupling is /i^""^ (/i^~'*~''ii 
for the two temperature model H), as can be seen from 
the dimensional analysis of sections O] and 4.2, this cou- 



pling is irrelevant in the renormalization group sense and 
therefore normally neglected. 

Under a Galilean transformation, to a reference frame 
moving with respect to the laboratory frame with velocity 
V, the coordinates and fields change according to 

r' = r-5ovt, (A.3) 

t'^t, (A.4) 

V'^,(r',t') - V'o(r-5ov^,t) =7Ao(r,t), (A.5) 

j[,(r', t') = j[,(r - go V t, t) - Jo(r, t) - v , (A.6) 

where the prime (') refers to parameters and dynamic vari- 
ables measured in the moving frame. 

Using these transformation laws, it is easy to show 
that the time and space derivatives of ■00 and Jq are trans- 
formed according to 



Vr' It' = Vr It , 

d . d . _ . 

^ Ir' = 7^ Ir +.go V • Vr It 



dt 



(A.7) 
(A.8) 



where |r etc. simply indicates which variable is being held 
constant when the derivative is taken. With these rela- 
tions, it is easy to show that the materia l der ivativ e d/d t 
which appears on the left-hand side of (A.l) and (A.2), 
i.e. 



d 
di 



dt ' 

d_ 

dt' 



-gojo • Vr It 

+30Jo • Vr' It' 



(A.9) 



is invariant under a Galilean transformation, i.e. it pre- 
serves its form on going fro m on e reference frame to an- 
other, as indicated in Eq. (A. 9). The right-hand side of 
Eqs. (A.l) and (A.2) can also be seen from Eqs. (A.5) to 
(A.7) to be trivially invariant, given the fact that v is a 
constant vector. It is thus shown that the Langevin equa- 
tions describing model H are invariant under a Galilean 
transformation (the distribution of the noise being the 
same in both reference frames). 

This invariance must be preserved under renormaliza- 
tion, i.e., when we substitute -^Oj Joi 50i • ■ -i by their renor- 
malized counterparts ip, ], g, etc. For this to happen, the 

and Z^ have to compensate 



renormalization factors Zg 

each other in Eq. ( lA.Qj ), i.e., one musu nave ^g-^j 
which is Eq. (3.14). Notice that the same reasoning also 



applies to the two-temperature model H, once go is sub- 



the left-hand side of (A_^). This term is normally not in- 



stituted by 50 = V D^/D^ g^ 
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